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We study time series concerning rare events. The occurrence of a 
rare event is depicted as a jump of constant intensity always occurring 
in the same direction, thereby generating an asymmetric diffusion pro- 
cess. We consider the case where the waiting time distribution is an 
inverse power law with index fi. We focus our attention on fx < 3, and 
we evaluate the scaling 5 of the resulting diffusion process. We prove 
that 5 gets its maximum value, 5=1, corresponding to the ballistic 
motion, at u = 2. We study the resulting diffusion process by means 
of joint use of the continuous time random walk and of the general- 
ized central limit theorem, as well as adopting a numerical treatment. 
We show that rendering asymmetric the diffusion process yields the 
significant benefit of enhancing the value of the scaling parameter 5. 
Furthermore, this scaling parameter becomes sensitive to the power 
index /x in the whole region 1 < u < 3. Finally, we show our method 
in action on real data concerning human heartbeat sequences. 
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Abstract 



1 Introduction 



One of the basic expectations of the statistical analysis of time series of bi- 
ological, sociological and financial interest is that the process under study 



might depart from fully random behavior. Following the prescriptions of the 
authors of Refs. 0, |3| we use the time series to create a sort of random walk 
trajectory whose statistical properties are then carefully studied to detect 
the deviations from ordinary Brownian motion. To do so a given site of the 
time series is interpreted as the "time" at which the walker makes a jump. 
The intensity of the jump is determined by the value of the time series at 
that given site, and the walker jumps either forward or backward according to 
whether the sign of the value is positive or negative. The most natural way to 
establish a comparison between this walker and the ordinary random walker 
of statistical mechanics would be to compare the two resulting diffusion pro- 
cesses to one another, but this would imply that infinitely identical walkers 
are available. Unfortunately a time series means only one trajectory, thereby 
generating the problem of how to create a number of "independent" trajec- 
tories large enough to produce a diffusion process. The authors of Refs. [Q, |3| 
addressed this problem with an ingenious technique called Detrended Fluc- 
tuation Analysis (DFA) which allows us to bypass this difficulty. The whole 
pseudo-random trajectory is divided into non-overlapping intervals of length 
I, each of which can be imagined as being a trajectory running for a "time" 
I. Thus we have a collection of a relatively large number of "independent" 
trajectories that can be used to produce the wanted diffusion process. These 
trajectories might be affected by a local trend and this local trend might 
result in an artificial deviation from ordinary Brownian motion. For this 
reasons the authors of Refs. P, [| for any interval evaluate the local bias and 
then record only the deviation from the local trend. The square root of the 
resulting variance, denoted by them with the symbol F(l), is expected to 
yield for very large Z's 

F(/)oc/ 5 , (1) 

where 5 denotes the scaling parameter, a crucial property to assess the de- 
parture from ordinary Brownian motion. In fact, ordinary Brownian motion 
should yield 5 b = 0.5 and the deviation of the observed value of 5 from 
5b = 0.5 is considered to be a measure of the deviation from a totally ran- 
dom behavior. 

It has to be pointed out that scaling is a property implying that 

pGM) = j^(£) (2) 

and that in some special cases, discussed in this paper, the scaling value 
detected by means of Eq.(P does not afford a satisfactory account of the 
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property of Eq.(^). Furthermore, as we shall see, the scaling 5, even if it 
reflects correctly the property of Eq.(fJ), might not disclose an exhaustive 
information on the dynamic source of the deviation from ordinary Brownian 
motion. 

Here we illustrate a technique of analysis, called Diffusion Entropy (DE) 
method, originally applied to detect memory effects in time series of sociologi- 
cal interest Q. We shall show that this technique can, in part, solve the earlier 
mentioned problems. We use a sequence of the same kind as that adopted 
in earlier work to mimic the statistical properties of DNA sequences ||. This 
sequence consists of strings of extended length, filled with only +'s or — 's. 
We consider as significant event only the change of sign. Thus, this kind of 
sequence is changed into a new sequence realized by assigning the vanishing 
value to the sites with a sign identical to that of the following one. Only the 
sites with a sign opposite to that of the next site are assigned a non vanish- 
ing value, which is always the same positive quantity, W. Thus we obtain 
a sequence of rare events. When one of these events occurs, the random 
walker makes a jump of constant intensity, always in the same direction. We 
prove that the resulting diffusion process makes much more ostensible the 
correlated nature of the original time series. 

We also show that this procedure makes faster the detection of the proper 
scaling parameter 5, namely, a scaling value which is a fair reflection of the 
property of Eq.(0). As earlier mentioned, the methods determining scaling 
through variances might yield scaling parameters significantly departing from 
the ideal scaling value. This is the scaling that an artificial time series, 
whose asymptotic behavior is dictated by compelling theoretical arguments, 
is expected to reach in the long-time limit. 

The outline of the paper is as follows. In Section ^| we describe the main 
properties of the DE method, and with the help of an artificial sequence 
we show that the observation of a symmetric diffusion process results in a 
slow convergence towards the emergence of the expected scaling parameter. 
Section [3] illustrates a theoretical approach to the scaling of the asymmetric 
diffusion process using two distinct procedures, the Continuous Time Ran- 
dom Walk (CTRW) and the Generalized Central Limit Theorem (GCLT), 
yielding the same result. This section also shows, with the help of the DE 
method applied to some artificial sequences, that the detection of the correct 
scaling is faster than in the case of symmetric diffusion. Finally, Section |] 
shows our detection technique at work on real data. The illustrative example 
of Section |] is based on a time series resulting from the record of human heart 
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beats. A balance on the results obtained in this paper is made in Section [5]. 



2 Diffusion Entropy 

The main idea of this approach to scaling is remarkably simple. Let us 
consider a sequence of M numbers , with i = 1, . . . , M. The purpose of 
the DE algorithm is to establish the possible existence of a scaling, either 
normal or anomalous, in the most efficient way as possible, without altering 
the data with any form of detrending. Let us select first of all an integer 
number I, fitting the condition 1 < I < M. This integer number will be 
referred to by us as "time". For any given time I we can find M — I + 1 
sub-sequences defined by 

^ s) =^+s, s = 0,...,M-l. (3) 

For any of these sub-sequences we build up a diffusion trajectory, labelled 
with the index s, defined by the position 

^(0 = Ee! s) = E^- ( 4 ) 

i=l i=l 

Let us imagine this position as referring to a Brownian particle that at 
regular intervals of time has been jumping forward or backward according 
to the prescription of the corresponding sub-sequence of Eq. (|3]) . This means 
that the particle before reaching the position that it holds at time I has been 
making I jumps. The jump made at the z-th step has the intensity | and 

(s) 

is forward or backward according to whether the number Q is positive or 
negative. 

We are now ready to evaluate the entropy of this diffusion process. To do 
that we have to partition the x-axis into cells of size e(Z). When this partition 
is made we have to label the cells. We count how many particles are found 
in the same cell at a given time /. We denote this number by Ni(l). Then we 
use this number to determine the probability that a particle can be found in 
the i-th cell at time Pi(l), by means of 



At this stage the entropy of the diffusion process at the time I is determined 
and reads 

S d {l) = -Y^PiWnM)]. (6) 

The easiest way to proceed with the choice of the cell size, e(l), is to assume 
it independent of I and determined by a suitable fraction of the square root 
of the variance of the fluctuation 

Before proceeding with the illustration of how the DE method works, it is 
worth making a comment on the way we use to define the trajectories. The 
method we are adopting is based on the idea of a moving window of size I that 
makes the s — th trajectory closely correlated to the next, the (s + 1) — th 
trajectory. The two trajectories have I — 1 values in common. The DFA, 
on the contrary, is based on non-overlapping windows, and consequently 
trajectories with different labels are totally independent the one from the 
other. A motivation for our choice is that we are having in mind a possible 
connection with the Kolmogorov Sinai (KS) entropy || [7fl. The KS entropy 
of a symbolic sequence is evaluated by moving a window of size I along 
the sequence. Any window position corresponds to a given combination of 
symbols, and, from the frequency of each combination, it is possible to derive 
the Shannon entropy S(l). The KS entropy is given by the asymptotic limit 
li l mi-+ 00 S{f)/1>- We believe that the same sequence, analyzed with the DE 
method, at the large values of / where S(l)/l approaches the KS value, must 
yield a well defined scaling 5. To realize this correspondence we carry out the 
determination of the DE using the same criterion of overlapping windows as 
that behind the KS entropy. 

Details on how to deal with the transition from the short-time regime, 
sensitive to the discrete nature of the process under study, to the long-time 
limit where both space an time can be perceived as continuous, are given 
in Ref. ||. Here we make the simplifying assumption of considering so large 
times as to make the continuous assumption valid. In this case the tra- 
jectories, built up with the above illustrated procedure, correspond to the 
following equation of motion: 

where denotes the value that the time series under study gets at the 
/ — th site. This means that the function £(Z) is depicted as a function of 
t, thought of as a continuous time t — I. In this case the Shannon entropy 



reads 

/oo 
dxp(x,t)ln\p(x,t)]. (8) 
-oo 

We assume 

This is a generalization of the ordinary scaling assumption of Eq.(^), which, 
in fact, can be recovered by assuming 5(t) to be time independent so that 
we can denote it as 5. For the sake of simplicity we make the simplifying 
assumption that F(y) maintains its form, namely, that the statistics of the 
process are unchanged. This is questionable in the case where the scaling is 
a function of time. However, making this assumption we obtain the following 
results. Let us plug Eq.(]^) into Eq.(||). After a simple algebra, we get: 

S(t)=A + 5(t)t, (10) 

/oo 
dyF(y)ln[F(y)} (11) 
-oo 



where 



and 

T = ln{t). (12) 

It is evident that this kind of technique to detect scaling does not imply 
any form of detrending, and this is one of the reasons why some attention 
should be devoted to it. An interesting way to check the efficiency of this 
technique is realized by the study of the artificial sequence of Ref. ||. This 
sequence is built up in such a way as to realize long sequences of either +'s or 
— 's. The probability of finding a sequence of only +'s or only — 's of length 
t is given by 

J7*-1 

^ = ^-VjrTry (13) 

Here we focus our attention on the condition \x < 3 and we raise the reader's 
attention on the interval [2,3]. In fact, this kind of sequence is the same as 
that adopted in earlier work [lOj for a dynamic derivation of Levy diffusion, 
which shows up when the condition 2 < fi < 3 applies. It corresponds to a 
particle travelling with constant velocity throughout the whole time interval 
corresponding to either only +'s or only -'s, and changing direction with no 
rest, at the end of any string with the same symbols. 
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We will refer to this model as Symmetric Velocity Model (SVM). We know 



from the theory of Ref. |10| that the scaling of the resulting diffusion process 
when 2 < \i < 3 is 

Note, however, that this diffusion process has a finite propagation front, with 
ballistic peaks showing up at both x = t and x = —t. The intensity of these 
peaks is proportional to the correlation function 

T~+i) ' (15) 

As a consequence of this fact, the whole distribution does not have a single 
rescaling. In fact, the distribution enclosed between the two peaks rescales 
with 5 of Eq.([L4]) while the peaks are associated to 5 = 1. Furthermore, it is 
well known |]10| that the scaling of the second moment is given by 



Thus, it is expected that the scaling detected by the DE method might 
not coincide with the prediction of Eq.(|TJ]) for the whole period of time 
corresponding to the presence of peaks of significant intensity. We think 
that the Levy scaling of Eq. (|TJP will show up at long times, when the peak 
intensity is significantly reduced. This conjecture seems to be supported by 
the numerical results illustrated in Fig.[l]. We see in fact that the scaling 
predicted by Eq.(THD is reached after an extended transient, of the order of 
about 20,000 in the scale of Fig.|l| This time interval is about 2000 larger 
than the value assigned to the parameter T, of Eq. flT3"|) , which is, in fact, in 
the case of Fig.|l|, T= 10. 

In conclusion, this section proves that the DE method applied to the SVM 
yields, for the scaling parameter 5, the correct value of Eq.(|14]), rather than 
the value that would be obtained measuring the variance of the diffusion 



process, Eq.(16). However, the time necessary to make this correct value 



emerge is very large. Furthermore, as proved by the theory of Ref. ||1 1|| , the 
adoption of SVM would make the scaling parameter 5 insensitive to [i in the 
whole interval [1,2]. This means that the adoption of the DE method would 
not allow us to distinguish a process with fi very close to 1 from one with n 
very close to 2. In Section ^| we shall see how to overcome these limitations. 
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3 Asymmetric Diffusion 



Let us consider again the kind of artificial sequence that we analyzed in 
Section §. Let us change perspective. Let us explicity focus on the events 
corresponding to changing sign. For this reason we build up a new sequence 
where the sites of the laminar region are given the vanishing value. This 
means that the states of constant velocity are here perceived as states of 
rest. The random walker can make a jump only at the moment when the 
time series analyzed in Section II changes signs. If at the moment of changing 
sign the walker makes a jump, of the same intensity, forward or backward, 
according to whether the sign change is positive or negative, we shall refer 
to this process as Symmetric Jump Model (SJM). The scaling corresponding 
to the SJM has been studied by Shlesinger in the pioneer paper of Ref.|12|, 



which yields the following prescription when the condition 1 < [i < 2 applies, 

5 = i0*-i). (17) 



When /i > 2, the theory of Ref. |12| predicts 



5 = 1/2. (18) 

We see that the whole region > 2, including that with \i < 3, would 
be indistinguishable from an ordinary Brownian diffusion. For this reason, 
in addition to focusing our attention on the events corresponding to the sign 
change, we decide to disregard whether the change of sign is from - to + 
or vice-versa. This generates a kind of asymmetric diffusion process that 
we term Asymmetric Jump Model (AJM). This section is devoted to the 
discussion of the AJM by means of two distinct theoretical approaches, both 
resulting in the same prediction, an interesting form of asymmetric Levy 
diffusion. Then we generate two artificial sequences corresponding to > 2 
and fi < 2, respectively, and we analyze them by means of the DE method. 

3.1 Continuous Time Random Walk 

The first theoretical approach rests on the Continuous Time Random Walk 



(CTRW) of Montroll and Weiss This means that the distance between 
two nearest neighbor events, Tj = ti + \ — U, is considered to be a random 
number, described by the probability density ip{r). Note that we are here 
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using the notation r rather than t to make clear to the reader the differ- 
ence between the absolute time t and the time distance between two nearest 
neighbor events. As earlier said, we assume that the random walker always 
moves in the same direction. Moreover, we also set the condition that the 
steps have constant length. Thus, the probability that a step of length x is 
taken is given by 

IL(x)=6(x-W). (19) 

Applying to this case the CTRW formalism of Montroll and Weiss Jjl^], we 
write for p(x,t), the probability of the walker being at position x at time t, 

p(x, t) = F^L^pfk, s) = F- 1 L- 1 \ ] —^ — ], (20) 

a l-A(Jfe)#0 

where X(k) and ip(s) are the Fourier and Laplace transform of U(x) and ip(t), 
respectively. Note that 

X(k) = exp(ikW). (21) 
By Taylor series expansion of Eq . (f2~0"|) we get: 



p(k, s) = - — ^ V ^ n e mkW (22) 



n=0 



By evaluating the inverse Fourier transform of Eq.(p2|) we arrive at: 

p(x,s) = ^^2^ n 5(x~nW). (23) 

S n=0 

This means that for large distances we can adopt the following expression 

p(x,8) = $(8)]*?—^. (24) 



We shed light on the meaning of Eq. fl2"3| ) with the following remarks. Let 
us denote by ip x (t) the probability of making x jumps up to time t, with the 
last jump occurring at t. Let us denote by $ n (i) the probability of making 
n jumps at earlier times, t' < t. The two functions are connected by the 
relation ^ 

* n (t) = [ Mt')*(t'-t)dt', (25) 
Jo 
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where 



/oo 
tP(t')dt'. (26) 



Using the properties x£) n (s) = [tjj(s)] n and $f(s) = (1 — ip(s))/s, and Laplace 
transforming Eq . (f25|) , we get 



KM = Ms)Hs) = [mr^—^-. (27) 



We see that the result of Eq.(p4|) can be recovered by identifying n with 
x/W, as it is legitimate to do when the number of jumps done is so large as 
to perceive n as a continuous number. 

Let us focus now our attention on the case where the distribution of 
waiting times, if)(t), has the following asymptotic behavior 

, , . const 

m » — , (28) 

with fi > 1 so as to insure the normalization condition. Let us focus our 
attention on the case where 

1< /I < 2. (29) 

This is the region where jumps of intensity x occurring with equal prob- 
ability in both the positive and the negative direction of the x-axis, i.e. the 
SJM, would result |TT| in a diffusion process with the scaling 6 of Eq.([l7|). 



We now show that, as a benefit of adopting the AJM rather than the SJM, 
we derive an enhanced scaling parameter. In fact, for very small values of s, 

tp(s) w 1 - cs 7 , (30) 

where 

7 = //-!. (31) 



This is the same property as that used by Zumofen and Klafter|IIJ. By 
plugging Eq.(^) into Eq.(^4j) and by inverting the Laplace transform with 
the nice method of ref.|p.4|, we obtain 



p(x,t)»-^£(4;>7,-7), (32) 

xi x~< 
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where L(y, 7, —7) denotes a fully asymmetric Levy stable law of index 7[I5|. 
It is straightforward to express this important result in the form of Eq.@ 
with the time independent scaling 

6 = ijl-1. (33) 

By comparing Eq. (|33| ) to Eq. ( |T7| ) it becomes evident the scaling enhancement 
produced by making the diffusion process asymmetric. In fact, at /x = 3/2 a 
transition is made from subdiffusion to superdiffusion. 

3.2 Generalized Central Limit Theorem 



The result of Section |3.1| can be recovered from within the attractive perspec- 
tive of the Generalized Central Limit Theorem (GCLT) ||16|| in a form derived 
from properly adopting the procedure of Feller JTTj to our purposes. To make 
closer the connection with the formalism of Feller let us assume W = 1. The 
function p(x, t) in this case is nothing but the probability for the walker to 
make x steps up to time t, or, equivalently, the probability of occurrence of 
x events up to time t. It has to be stressed that according to the perspective 
expressed by Refs. jl8|, [0| , in accordance with the arguments of Gaspard and 
WangPP[|, this way of proceeding can also be used to measure the complex- 



ity or randomness of the process under study, in the special case where the 
time series under study is totally memory- less. In fact, in this specific case 
the occurrence of any sporadic event does not have any connection with the 
earlier or subsequent events. It has to be pointed out that the kind of time 
series that we are having in mind in this paper is characterized by events of 
the same intensity, or made so by ignoring the real intensity of each event. 
If there is no correlation between intensity and time occurrence of the event, 
the simplifying assumption that the events have the same intensity does not 
have any significant consequence on the scaling of the resulting diffusion. 

If the simplifying assumption applies, that no form of memory exists, in 
addition to the inverse power law nature of the waiting time distribution ip(t), 



we can safely apply the Feller arguments. Feller [171 studied the case where 



the survival probability ^(t) has the following time asymptotic expression 

*(*) « ^ (34) 
The Feller notation is related to our notation by 7 = \i — 1. In the article 



of Ref.[17|, Feller studies bot the case < 7 < 1 and the case 1 < 7 < 2, 
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which refers to our two regimes 1 < /x < 2 and 2 < 7 < 3, respectively. In 
the former case, Feller finds the important property 

Pr (N t > wf) = G 7 (w~ 1/7 ), (35) 

where, according to the notation used by Feller |J7] Pr(N t > y) denotes the 
probability that the number of events occurring at times t' < t is greater than 
y. The function G 7 (x) denotes the stable law whose characteristics function 
is given by 

G y (z) = exp |-|2| 7 cos (-yj - % sin (J^-) sgn(z) r(l-7)|. (36) 

The symbol T denotes the Euler function. The meaning of the function 
G 7 (w) is made transparent by the case where the walker makes jumps in 
the same direction, with the probability U(x) oc 1/x 7 , at regular interval of 
times. In this specific case, which is the case of an asymmetric Levy flight, 
the position of the random walker at time t, X t , fulfills the property 

Pr (X t < wt lh ) = G 7 (w). (37) 



The function on the right hand side of Eq. ( |37j ) is the same as that appearing 
in the right hand side of Eq. (|35|). However, we know that it has to do with 
an asymmetric Levy process. Consequently, it is a stable law that the theory 
of [15|] makes it possible to write under the form 



G 7 H = J L(y,-y,-j)dy, (38) 


where L(y, 7, —7) is the function defined in Section |3.1|. We note that Eq. 



35p can also be written as 

-00 

dN t p(N t ) = G 7 (w 7 ). (39) 



By differentiating Eq. ( |39D with respect to w and settingA^ = x (the position 
of the walker at time £, labeled clS OC ^ IS equal to the number of steps taken 
up to time t, namely, N t ) we get 

p(x,t)« * £(^,7,-7), (40) 

x 1 X"i 
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which agrees with the scaling predicted by the CTRW arguments used in 



section 3.1 since it leads to the scaling parameter 5 = fi — 1. 



In the case where 1 < 7 < 2 we are led by the corresponding Feller pre- 



scription |L7| to express the probability distribution in the reference trans- 
lating with the mean velocity (the inverse of the mean waiting time). With 
an argument similar to the one earlier used we obtain 

p(x,t) » — L(4>7,-7), (41) 
ti t~t 

and it is now easy to prove that this yields, for 2 < /x < 3, the scaling 
prescription 

, = i^T- < 42 > 

Results are summarized in Tab. 1. 



3.3 Numerical Calculations 

This section is devoted to illustrating the results of numerical calculations 
that confirm the benefits of the strategy suggested in this paper. First of all 
with Fig.^ we show that the AJM makes it possible to reach much earlier 
the correct Levy scaling of Eq. (|I3D . The numerical work shows that the time 
scale at which the DE method yields the correct scaling is, in this case, equal 
to about 2,000 (T = 10). 

We also prove that the resulting scaling value is a genuine reflection of 
the property of Eq.@. To do so, we express the diffusion process in the 
reference framework moving with the velocity TW/ (fi — 2) corresponding to 
the theoretical prediction of section p.2| . Then we consider a reference time t no 
and later times t n > t no . For any time t n we consider the squeezing parameter 
R(n) = (t n /t no ) s , with 5 given by the scaling detected by the DE method. 
The distance scale is squeezed by the quantity R(n) and the ordinate scale 
is enhanced by the quantity 1/R{n). If the scaling detected by the DE is 
correct and if this procedure is applied to the distribution corresponding to 
the time t n this scaling procedure should make it identical to the distribution 
corresponding to the time t no . Fig. |^ shows that this is really so, since all the 
distributions at times t n are made to virtually coincide with the distribution 
at time t no . 

We have explored also the case /i < 2. Figure f| shows the result of the 
DE method in this case, and Fig.|5] shows the result of scaling back to earlier 
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times the distributions corresponding to later times, in this case. The result 
is not as accurate as that illustrated by Fig|| However, we judge it to be 
satisfactory, especially if we take into account the fact that the region /i < 2 
corresponds to the case of extremely large time correlations, with an infinite 
sojourn time. 



4 The Method in Action in the Case of Real 
Sequences 

As an example of how this technique can be applied to real data we con- 
sidered time series concerning human heartbeats. These series were taken 
from Physionet EI] and are a record of time intervals between successive 
beats. To apply our method we adopted a coarse graining procedure: times 
between beats were first multiplied by 10, with only the integer part of the 
resulting number taken into account. This sequence of integer numbers was 
then transformed into a new one according to the following prescription: we 
wrote if the number was equal to the preceeding one, 1 if it was differ- 
ent. We note that this would have resulted almost in a sequence of only l's, 
without any 0, had we not adopted the coarse graining procedure. From this 
sequence we generated an asymmetric diffusion process (AJM) which we then 
studied by means of DE method described in the previous sections to derive 
the scaling parameter 5. This corresponds to having a walker taking a step in 
a fixed direction when the time interval between successive beats, as seen on 
a coarse grained scale, differs from the preceding one. We also analyzed the 
probability distribution tp(t) of waiting times between the l's, which shows 
an inverse power-law behavior with index /j > 2. We thus expect the correct 
scaling to be predicted by Eq.(^) rather than by Eq. fl33|) . Results obtained 
from different databases are reported in table |5|, whose third column repre- 
sents the value of /x we derive from 5. An example for a particular database 
(Ltdbl5814) is shown in Fig.|5| where we derive the scaling 5 using DE, and 
in Fig.^, where we show how the inverse power law of index /i = 1 + 1/5 
indeed fits the data used to obtain if)(t). 
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5 Concluding Remarks 



This paper is a sequel to that of Ref . || , where the method of DE was applied 
for the first time. However, the emphasis here is on the detection of the 
scaling more than on the property that the authors of Ref.[|J referred to 
as memory. On the basis of the results of the present paper, the memory of 
Ref. [0 seems to have to do with the transition from the short-time to the long- 
time regime. The main result here reached, however, has to do with proving 
the convenience of recording the times at which events occur, an event being 
a sign change, regardless of whether it implies an increase or a decrease of 
the intensity signal. If we record the time at which the event occurs, and 
we interpret it as a jump of the random walker always in the same direction, 
with a fixed intensity, then the diffusion process becomes asymmetric, with 
the big benefit of enhancing the value of the scaling parameter in the interval 
1 < /i < 2, and of making it sensitive to the power law index /x in the whole 
interval 1 < fi < 3. The AJM has to be compared to the SJM, a model that 
the theoretical analysis of the 1974 pioneer work of Shlesinger [12| shows to be 
insensitive to /i in the range 2 < fi < 3. The value of the scaling parameter 5, 
stemming from the AJM, as earlier repeatedly pointed out, is enhanced with 
respect to that of the SJM, in the whole interval 2 < fi < 3. If we compare 



the AJM to the SVM of Ref.|TTJ, we find that it provides information on /i 
even in the interval 1 < /i < 2, where the SVM would steadily result in the 
ballistic prescription 5 = 1. 

It is interesting to notice that the AJM proposed by this paper emphasizes 
the phase transition character of a transformation moving from fi > 2 to 



fi < 2. As pointed out in Ref. ||19||, this is a transition from a stationary to a 



non- stationary regime. In fact, Ref. |B| shows that /i < 2 implies the lack 
of an invariant distribution, thereby making any observation carried out for 
a finite time equivalent to monitoring an out of equilibrium process. We are 
convinced that the systematic application of the prescriptions of this paper 
to studying time series generated by complex dynamic processes, might lead 
to discover these non-stationary cases. 

The application to the real sequences of physiological interest of Section 
[| shows that the detection of the scaling parameter 5 is much more accurate 
than the direct evaluation of \i from the experimental data, the only ambigu- 
ity being the fact that the value of 5 is compatible with both /i(l) = 1 + 1/5 
and /i(2) = 1 + 5. So, the experimental error affecting the direct determi- 
nation of delta must only be smaller than the difference /i(l) — /u(2). Only 
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when the value of fi is very close to the transition value fi = 2, the adoption 
of the DE method to detect fi becomes ambiguous, and additional criteria 
must be found. 

A possible efficient criterion could be given by the joint use of the SVM 
and AJM method. In fact, as clearly illustrated by Table [l], if the SVM 
yields 5 = 1, this means that fi < 2. If, on the contrary, the SVM, within 
the statistical accuracy of the numerical procedure, yields the same scaling 
as the SVM, we have to select ji > 2. 

Another benefit of the adoption of the AJM is the faster attainment of the 
correct scaling. It is not yet quite clear which is the true reason why the cor- 
rect scaling of (0), in the region 2 < fi < 3, is significantly faster than with 
the SVM. One reason might be that the probability distribution of the AJM 
diffusion process is sharper than that of the SVM. Some more research work 
is required to assess this point. We hope that the results already obtained 
are interesting enough as to draw the attention of the researches of this field 
of investigation on the benefit stemming from converting the experimental 
data into asymmetric diffusion processes. As the last but not least remark, 
we want to point out that the DE method is probably the most convenient 
way to detect the scaling of an asymmetric diffusion processes, since this 
technique is not affected by the presence of biases, and does not require the 
adoption of any form of de-trending. 
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Table 1: The scaling parameter 5 as a function of for SVM, SJM and 
AJM. The AJM is the only model where 5 always changes upon changing \x. 



Archive 


S 


M= 1 + 1/5 


14046 


0.777 ±0.002 


2.29 


14134 


0.863 ±0.006 


2.16 


14172 


0.876 ±0.002 


2.14 


14184 


0.804 ±0.005 


2.24 


14157 


0.855 ±0.010 


2.17 


15814 


0.872 ±0.002 


2.15 



Table 2: Values of the scaling parameter 6 obtained with the DE method. 
The archive numbers correspond to the denomination of MIT-BIH Long- 



Term ECG Database as reported in ref . [BT1 . 
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Figure 1: The diffusion entropy as a function of time. The numerical 
method is applied to the artificial sequence of Section ||, with /i = 2.5, studied 
according to the SVM prescription. According to the theoretical arguments 
of the text the scaling parameter 5 is the slope of the straight line fitting the 
numerical results at large times, which yields in this case 6 = 2/3 = — 1) 
(T = 10). 
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Figure 2: The diffusion entropy as a function of time. The numerical 
method is applied to the artificial sequence of Section [|, with \i = 2.5, studied 
according to the AJM prescription. According to the theoretical arguments 
of the text the scaling parameter 5 is the slope of the straight line fitting the 
numerical results at large times, which yields in this case 5 = 2/3 = l/(/i — 1) 
(T = 10). 
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Figure 3: Scaled distributions as a function of x. These scaled distributions 
refer to the artificial sequence of Section 0, with \i = 2.5, dealt with according 
to the AJM prescription. The times considered are: t n =2500 '+', 2000 V, 
1500 'x', 1000 Note that t m = 1000. Time is expressed in units of T, 
with T defined by Eq.(|13|), and T = 10. 
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Figure 4: The diffusion entropy as a function of time. The numerical 
method is applied to the artificial sequence of Section 0, with \i = 1.8, studied 
according to the AJM prescription. According to the theoretical arguments 
of the text the scaling parameter 5 is the slope of the straight line fitting the 
numerical results at large times, which yields in this case 5 = // — 1 = 0.8 
(T = 10). 
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Figure 5: Scaled distributions as a function of x. These scaled distributions 
refer to the artificial sequence of Section with \i = 1.8, dealt with according 
to the AJM prescription. The times considered are: t n = 2500 '+', 2000 V, 
1500 'x', 1000 Note that t no = 1000. Time is expressed in units of T, 
with T defined by Eq. (|IJ), and T = 10. 
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Figure 6: The DE method applied to the heart beating sequence Ltdbl5814, 
changed into a diffusion process according to the AJM prescription. The 
straight line is the best fit made with Eq. QIQp and it yields 5 = 0.872. 
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Figure 7: Coarse-grained waiting time distribution ip(t) for Ltdbl5814. The 
straight line is an inverse powe law curve of the form C/f 1 . According to the 
theory of section j|, \i = 1 + 1/8. 
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